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LONG TIME CONVERGENCE FOR PHASE FIELD MODELS 



1 Introduction 

The present note is devoted to the analysis of the regularity and long-time behavior 
properties of the following class of PDE's systems modeling phase change phenomena 



Here, the unknowns are the relative temperature d (i.e., some critical freezing or melting 
temperature r c has been normalized to 0) and the order parameter, or phase field, X; 
both are functions of the spatial variable x (ranging in a bounded, connected, and 
sufficiently smooth domain Q C R d , 1 < d < 3) and of the time t e [0, oo) (let us use 
the shorter notation oo in place of +oo). In the above system, e, 5 > are relaxation 
parameters; the function W : domH 7 — > M represents a configuration potential in X 
and is assumed to be convex in its principal part; A(-) is a possibly non linear function 
with ed + X(X) yielding the internal energy of the system; finally, / stands for the 
heat source. However, the main novelty in (ll.ip - (ll.2l) is given by the presence of the 
nonlinear, but convex function j. This nonlinearity can be justified both on the physical 
and on the analytical side. Actually, several well-known models can be included in this 
general framework. For example, j(r) = r 2 /2, corresponds to the Caginalp system 
[6], while j(r) = — log(r + r c ) + r/r c gives the so-called Penrose-Fife model [2U [25] 
(the complication of the expression, which is not the usual one, is compensated by the 
nice property j'(0) = 0). Furthermore, also intermediate choices for j correspond to 
meaningful cases: for instance, a combination of the previous two expressions provides 
a variant of the Penrose-Fife model with special heat flux law introduced in [TOl [11] . 

Various systems concerned with the abovementioned phase field models, possibly 
including non-smooth potentials and also applying to martensitic phase transforma- 
tions, have been intensively investigated in the last years. Among a number of recent 
contributions, let us quote and [E] which treat the case of Neumann boundary con- 
ditions for the temperature; [22] devising a general (convex) framework for the study 
of Penrose-Fife systems; [TJ] involved with the analysis of the quasistationary (i.e., 
5 = in (11.21) ) Penrose-Fife model; [Tjl] that deals with the long-time behaviour and 
the study of inertial sets; [1] and [2] studying the long-time convergence of the Caginalp 
model with and without memory effects, and [30, EH] addressing the same questions 
for the Penrose-Fife model; [261 EZ] showing the existence of a uniform attractor for 
Penrose-Fife systems; [9] and [2H] in which a hyperbolic dynamics for X, characterized 
by an extra inertial term pX u in the left hand side of equation (11.21) . is considered. 

In all this framework, the occurrence of a general convex function j seems to be 
an interesting isssue, worth to be deepened. Recently, well-posedness of an initial and 
boundary value problem for (ll.ip - (ll.2l) has been shown in [29J for completely arbitrary 
(convex) j. More precisely, in [29] the homogeneous Neumann problem related to (II. ID - 
OL 21) is studied for a linear function A. We remark that, while the Neumann boundary 
conditions for the phase field X appear to be the most natural choice for phase field 
models, the case of no-flux conditions for was motivated in [29] by the purpose 
of studying some singular limits of the system. In particolar, in |29j it is shown that 
(ll.ll) - (ll.2p gives rise to the Cahn-Hilliard equation in the viscous form if e is sent to 



e0 t + \(X) t -Aj\#) = f, 
5X t - AX + W'(X) = X'(X) 



(1.1) 
(1.2) 
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and in the standard form if both e, 5 go to 0. On the other hand, the choice of the 
no-flux conditions for gives rise to some additional difficulties in the analysis that 
forced the author of [29J to restrict the range of the admissible potentials W in order 
to get some a priori control of the spatial average of the unknowns (see [291 assumption 
(28)], see also [8]). 

In this paper, we first extend the well-posedness result of [2U| . by adapting it to 
our slightly different setting. Then, we investigate here some further properties of the 
solution to suitable initial-boundary value problems related to fll.ip - fll.2l) . Namely, 
we shall concentrate our attention on the long-time behavior of the system from the 
point of view of cu-limits of solution trajectories. Since we are not interested in singular 
limits, we shall take e = 5 = 1 in the sequel. 

In our approach, the basic observation is that the system, at least in the case of no 
external source, admits the Liapounov functional 

X) := J (||VX| 2 + W(X)+j(#j), (1.3) 

which is obtained testing f ll.ip by j'l'd) and ( 1 1.2ft by X t , then taking the sum. Hence, if 
we introduce the new variable e := d + A(X) (internal energy), f ll.ip -( jl~2l) can be seen 
(at least in the case of no-flux conditions and no external source) as a gradient flow 
problem 

et = -d*, e £, X t = -d Ha £. (1.4) 

Here, the symbol d denotes here (sub) differentiation, and, more precisely, dn, x stands 
for the sub differential w.r.t. X in H := L 2 (Q) while <9* je indicates differentiation in e in 
the space if x (f2)* (which gives rise to the Laplacian in fll.ip ). 

In view of this variational structure, it is reasonable to expect good asymptotic 
properties of the solution as t goes to infinity. This has already been noticed in some 
particular cases. For instance, for the Caginalp model corresponding to j(r) = r 2 /2, 
the long time analysis has been performed in j2j [16] for various types of boundary con- 
ditions and assumptions on A, W. Instead, the Penrose-Fife case (j(r) = — logr) has 
been studied in [14] . referring only to the case of (nonhomogeneous) Dirichlet bound- 
ary conditions for the temperature. Concerning [T4], we point out that the problem 
addressed there presents a number of mathematical difficulties due the character of 
j(-) = — log(-), which is both singular at and unbounded from below. Actually, the 
precompactness of trajectories, leading to the existence of a nonempty cu-limit set, has 
been proved in p3] by strongly relying on ad hoc techniques to overcome, in particular, 
the noncoercive character of j at +oo. 

Based on these considerations, it seems very difficult to perform a long time analysis 
of f ll.ll) -( fl~2l) by taking j completely general as in [29]. Rather, we find there is room 



for a bit of compromise and ask at least that j has some coercivity: see ( hp j 1 ) below. 
On the other hand, j can exhibit a singular character too. Still for the purpose of 
coercivity, we also need boundary conditions for $ which give some (uniform in time) 
control on its space average (differently from the no-flux case of [29J). Then, we shall 
consider two sets of boundary conditions: 

i) $ = $oo, with /(tfoo) = 0, and<9 n X = 0, (1.5) 
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which we refer to as Dirichlet boundary conditions, and 

ii) - d n j\$) = V (j'(#) - f(&r)) and d n X = 0, (1.6) 

which we refer to as Robin boundary conditions, where represents the extremal 
boundary temperature and r\ is a positive constant. We point out that the latter case 
subsumes the presence of some (boundary) source term. We will see that this does not 
destroy the variational character of the system, nor does this the presence of a nonzero 
volumetric heat source / in ( 11. II) . provided / and the boundary datum are globally L 2 
in time. In our analysis, we will be able to consider nonlinear latent heat functions A 
and, what is more important, rather general potentials W with the only restriction that 
they should not exhibit minima at the boundary of their domain (like instead it could 
happen in the case of a double-obstacle potential). These could neither be considered 
in [H] , essentially for technical reasons, nor in [29] due to the quoted difficulties coming 
from the no-flux conditions for j'ty). We stress in particular that singular potentials 
W (i.e. which are +oo outside an interval I C R) are not completely easy to address 
even in the simpler case of the Caginalp model [TB] . 

In the framework of the long time analysis, we prove the following results. First, 
we demonstrate that the cu-limit of any solution trajectory is not empty and consists 
only of steady-state solutions ($oo,^oo)- More in detail, the component coincides 
with the (constant) temperature on the boundary (i.e., with the unique value such that 



j'($<x>) = 0, cf. assumption ( hpj'2 ) below); so that we also obtain the limit $(t) — > 
holds in a suitable sense as t /* 00 and not only for a subsequence. Conversely, in 
general we only have that, as t /* 00, X(t) is precompact in a suitable topology and 
any of its limit points is a steady state solution of ( 11.21) . In fact, ( 11.21) may have infinitely 
many stationary states due to nonconvexity of W. We also point out that a careful 
use of parabolic regularization effects is a key step in the proof of the convergence 
result. Our second theorem, which follows the lines of some recent work devoted to 
the large time behavior of degenerate parabolic equations P2] and phase transition 
systems (see, among others, [21 CE31 EH]), gives a sufficient condition under which the 
tu-limit consists of only one point. Namely, we can prove that, in the case when W 
is analytic in the subdomain where the solution component X lives, then the tu-limit 
set is a singleton and consequently also the entire trajectory of X(t) converges to X^. 
Here, the main ingredient of the proof is the so-called Simon-Lojasiewicz inequality 
[2"31 I3T] . that was originally [23] stated as a nontrivial (local) growth estimate for 
analytic functions of several complex variables in the neighbourhood of a critical point: 
Simon contributed by extending this inequality to the infinite dimensional setting, thus 
allowing to characterize the large time behavior of evolution systems with analytic 
nonlinearities. The last result we present is a convergence result which also establishes 
a (say, polynomial) rate of convergence for the L 2 (f2)-norm of X(t) — X^. In proving 
such a result, we follow a method from [7] and argue partly as in [T7j . 

Here is the plan of the paper. In Section [2] we introduce the functional framework, 
state precise hypotheses on the data and formulate our main results. Section [3] is 
devoted to the proof of the well-posedness of the system (Theorem 12.21) and of the 
basic uniform estimates (Theorem I2.6p . Finally, Section H] is concerned with all the 
properties of the cj-limit set (proof of Theorem 12.71 and Theorem I2.13P . 
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2 Main results 



Let Q be a C 1 ' 1 , bounded, and connected domain in M. d , 1 < d < 3, and let T := dQ. 
Set i7 := L 2 (Q), V := if 1 (£7), both endowed with their standard scalar products and 
norms. The norms in H and in H d will be equally indicated by | • | and (•, •) will denote 
the corresponding scalar products. Let also V := Hq(CI), endowed with the norm 
II ' \\v := |V • |. The symbol || • \\x will stand for the norm in the generic Banach space 
X and x*( m , -)x will denote the duality between X and its topological dual space X*. 
The space H will be identified with its dual, so that we have the chains of continuous 
embeddings V C H C V* and V C H C Vq. We introduce the elliptic operator 



A : V -> V\ 



v *(Av,z)v ■-- 



Vf • V-z, for v, z G V. 



(2.1) 



We also let J, J be open intervals of M, with 6/, and let I be an open and bounded 
interval containing and whose closure is contained in /. Then, our basic hypotheses 
on the data are 



A G CiJ(M), 3 A > : |A"(r)| < A for a.e. r6l, 

W G CK(J; [0, +oo)), 3 k > : W"{r) > -k for a.e. re I, 

3/i > : W'{r)/r > /x for a.e. r e I\I , 

j G CjJ(J; [0, +oo)), 3 a > : j"(r) > tr for a.e. r G J, 

3$oo G J : j"(^oo)=0. 



(hpA) 

(hpWl) 
(hp^2) 

(hpjl) 
(hpj'2) 



Let us spend some words on the hypotheses on W and j. Formula Qhpiyip says that W 
can be nonconvex, but just up to a quadratic perturbation (actually, r i— > W(r) + Kr 2 /2 
is convex). Instead, j is uniformly strictly convex by ( hp j 1 ) ; moreover, it has some 



coercivity property in the sense that, by ( hpj2 ), it attains its minimum value (which 
can be fixed at 0, for simplicity) at the point G J. The role of ( hpVT2 ) will be 
outlined later on. 

In the sequel, both W and j will be extended to the whole real line by means of the 
following procedure. First, we prolong j (resp. W + kH 2 /2) giving it the value +00 
outside J (resp. J); then, we take the lower semicontinuous regularization, and, finally, 
only from W we subtract nld 2 /2. This means, in the simpler (because convex) case 
of j that, if J is bounded and j "explodes" (taking the value +00) at its boundary, we 
simply extend it at +00 outside J. If, instead, J is bounded but j does not explode at 
least on one side of J, then we first close the graph of j and then extend it at +00. For 
W the procedure is slightly more complicated due to its possibly non convex character. 
In any case, the extended j and W + kH 2 /2 are convex and lower semicontinuous 
functions from R to [0, +00], so that their IR-subdifferentials are maximal monotone 
graphs coinciding, respectively on J and on J, with the "original" functions j' and 
W' + nld. 

With all these conventions in mind, our assumptions on the initial data are 



$oeh, j(tfo) e L\n), 
x G v, w(x ) g l 1 ^). 



(hp^o) 
(hpXo) 
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Remark 2.1. Recalling fll .3j) . let us notice that ( hp^ )-( hpX ) are equivalent to asking 
£(t9cb^o) < +00, i.e., that the initial energy is finite. Indeed, the property $ E H 



follows from the quadratic growth of j (cf. ( hp j 1 ) ) . Note also that, by (hpWl ) and 



Qhpj'ip , W and j assume only non-negative values and consequently the functional £ 
is non-negative. 



Concerning the heat source, we assume in any case 

fEL 2 (0,oc;H), 

but this hypothesis will need some refinement in the sequel. 
In case we work with the Robin boundary condition 

- aj'{0) = ri{j'(#) - j'(tfr)) and d n X = 0, 



(hp/) 



(2.2) 



where n is some positive constant, we also suppose that the external boundary tem- 
perature satisfies 



$r : r x (0,oo) 



measurable, j'($r) G L 2 (0, 00; L 2 (T)). 



(hp^i 



We remark that, by ( hp j 1 ) (which ensures the Lipschitz continuity of (j') x ) and 



( hpj2 ), (hp#r) entails that, for n / 00 (n denoting here the time variable), 

j'(Mn + -)) ->0, ^ r (n + -)^^oo strongly in L 2 (0, 1; L 2 (r)). (2.3) 
Next, we introduce the operator 



R : V -> V* 



v*(Rv,z) v = I Vt> ■ Vz + I vz. 

In Jr 



(2.4) 



Of course, R turns out to be the Riesz isomorphism associated with the (equivalent) 
norm on V given by 



\°\\r — v*(Rv,v) v = I I V-^l^ + 77 / v £ 

>n Jr 



(2.5) 



At this point, we are able to state our first result, related to existence and uniqueness 
of solutions (with e — 1), which slightly extends [29l Thm. 1]: 



Theorem 2.2. Let dhpAj) , ^pW^ - ^pW^ , dhpiy , dhpXj ), ( |hp7| ) 

nold. Moreover, assume either (Dirichlet conditions) 



B := -A : V — > r o *, V := Vb, 



IVoj 



:= /, 



(Dirichlet) 



or (JhpWrJ) and (Robin conditions) 



B :=R, V:=V, \\ ■ \\ v := \\ ■ \\ R , v *{g,v) v := (f,v) + V j'($ r )v, (Robin) 
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the last relation holding for allvEV = V, a.e. in (0, oo). Then, there exist a constant 

c > 0, depending only on the data A, W, j, $o, X , and a unique pair X) such that 
ifu = f(#) 

II^Hl 2 (0,oo;V*) + ||0|U»(O,oa;.H) < C, (2.6) 

||m||l2(0,oo;V) < C, (2.7) 

||^-^oo||l 2 (0,oo;V) <C, (2.8) 

IWU 2 (0,oo;ff) + ||X|| L oo( 0)OO; y) < C, (2.9) 

X £ C°([0, t]; V) n L 2 (0, t; # 2 (ft)) for all t > 0, and tf, X, u satisfy 

tf t + \(X) t + Bu = g inV*, (2.10) 

X t + AX + W'(X) = X'(X)u inV*, (2.11) 

a.e. in (0, oo), as well as 

0|t=o = 0o, X|t=o = X o- (2.12) 

Remark 2.3. We point out that, in the (IRobinl) case, ( |hp/D and ( hp-^p ) entail 

^GL 2 (0,oo;V*). (2.13) 

This is of course true also in the ( IDirichletl) case, in which we even have by assumption 
the better relation g £ L 2 (0, oo; H). 

We shall not give the full proof of Theorem 12.21 since it is just a variant of the proof 
given in [29|. More precisely, the same argument of [221 Sec. 3] can be used to obtain 
existence of a solution to a suitable regularization of the problem. Concerning the 
a priori estimates which are required to remove the approximation, some points are 
technically different from [291 Sec. 3] especially due to our choice(s) of boundary condi- 
tions. Nevertheless, the required estimates will be easily obtainable from the uniform 
bounds we shall prove in the long-time analysis. Finally, the compactness argument 
necessary to pass to the limit and the proof of uniqueness will be briefly sketched in 
the next section. 

Let us come now to our statement on further regularity properties of solutions. To 
this aim, we prepare an auxiliary result. 

Lemma 2.4. Let X be a Banach space and let 7 £ L 2 (0, oo;X) satisfy, for some 
V e [l,oo], 

7 £ W lj, (0,t;X) Vt£(0,oo), sup H\\ L p(t,t+i;X) < 00. (2.14) 

te[o,oo) 

Then, it turns out that 7 £ L°°(0, 00; X). Moreover, if p > 1, the strong convergence 
7 (i) — > holds in X , as t /* 00. 

Proof. First of all, the L°° bound can be shown using the Fundamental Theorem of 
Calculus, by a simple contradiction argument. Concerning the convergence to 0, let 
p > 1 and {t n } be an arbitrary diverging sequence of times. Setting 

7n(*) :=7(*» + *), te[0,l], 
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it is clear that 7n — > strongly in L 2 (0, 1;X). Thus, there is a subsequence, not 
relabelled, such that 7n (t) — > strongly in X for a.e. t G (0, 1). In particular, we can 
find a sequence {5k} C (0, 1), with 5k \ and such that, in X, 

lim j n (5 k ) = VkeN. 

n— >oo 

We then deduce from (12.141) 

h(tn)\\x < \\l(t n ) - 7 (t n + 5 k )\\x + \\l(t n + < C^ /P * + ||7n(4)||x, 

where p* is the conjugate exponent to p, and the latter quantity can be made arbitrarily 
small for n large enough. Due to the arbitrariness of {t n }, this shows that 7 (t) — > in 
X as t oo. | 

Remark 2.5. We point out that the second part of the statement above is false (even 
for X = M.) if one takes p — 1. Indeed, set, for n > 2, v n : R — > R given by 

v n (t) := n 2 X(_i/ n 2^ — n 2 A'( 0il / n 2), 

denoting the characteristic function, and define 



lit) : =^2 v n(t-n), t 



n=2 



g(t) := / 7 (s) ds 



Then, it is clear that # stays in L 2 (0, oo) nL°°(0, oo) and satisfies fl2~T4l) (with X = R). 
However, g(t) does not tend to for t /* oo. 

Theorem 2.6. Let flhpg) (|hipW TJ) ( |hpV^2[), JEpjTl )-( |Epj2| ), dEp^oT ), (JhpX^), JEpJI l 
and either ( IDirichletl) . or ( hp^r ) and ( IRobinl) . nold. Let also (12.141) hold for 7 = g, 
X = V* and some p G [l,oo]. Then, for all s > there exists a constant c > 0, 
depending on A, IV, j, $ ? ^0 an d s , sucn that 



SUp 11^11x2(^+1.^) + ||i?||l°o( S)00 . v) < C, 
i>s 

IMU°°(s,oo;V) < C, With It = /($), 

||^t|U oo (a,oo;f0 + SU P ll"^i|U 2 (i,m;V) + ||^IU°°(s,oo;iT 2 (r2)) < C, 
t>s 

\\W'(X)\\ Loai8t0O . H ) <c. 
If, additionally, 

g t eL«(0,oo;V*) 
for some q G [1,2], then we also have 

\Wt\\L 2 {s,oo;H) + ||^|U 2 (s,oo;V) < C. 



(2.15) 

(2.16) 
(2.17) 

(2.18) 
(2.19) 

(2.20) 



In particular, the above theorem provides a priori estimates which have a uniform 
character for large times. Of course, this uniformity cannot be proved for the L 2 in 
time norm of X in ( 12.9ft . 

Let us now move to the study of long-time behavior, starting from existence of a 
nonempty cu-limit set. 
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Theorem 2.7. Let ( [hpA] , flhpWlD - flEpIFSl , JEpjU - jEpj^D , dhp^j ), flhpXj ), flEpJl ) 
and either ( IDirichlet| . or ( hp-$ r ) and (IRobinP . nold. Let aiso ( 12.141) hold for 7 = g, 
X = V* and some p G [1, 00]. Moreover, let us assume that, either 



g t eL 2 (0,oo;V*) 



or 



3c>0, a<3: \j"(r)\ < c(l + \j'{r)\ a ) Vr G J. 
Then, we have that, as t /* 00, 

u(t) — >■ 0, #(t) — > weakly in V and strongly in H. 



(2.21) 



(2.22) 



(2.23) 



Moreover, any diverging sequence {t n } C (0, cxo) admits a subsequence, not relabelled, 
such that 

X(t n ) -> Xoo weakly in H 2 (Vt) and strongly in V n C '(H) , (2.24) 
wnere Xoo is a solution of the stationary problem 

AXoo + W(Xoo) = in V*. (2.25) 



Remark 2.8. Assumption (12.141) on (7 of course reinforces the convergence properties 
of g to (cf. Qhp/P and, in the flRobinl) case, (hp'flrD)- 



Remark 2.9. The property (12.221) is actually not very strong. For instance, in the 
situation when J = R, then ( 12.22ft is fulfilled provided j has a polynomial, or even 
exponential, growth at infinity. Moreover, ( 12. 22ft also holds for the laws mentioned 
in the Introduction (corresponding to combinations of the Caginalp and Penrose-Fife 



models). A case in which ( I2.22p does not hold (but ( hpj'l )— ( hpj'2 ) do hold) is given by 



j'(r) = r + log(l + r/r c ), due to the singular behavior of j' in proximity of — r c . 

Our last result characterizes the u;-limit as a singleton in case the potential W is 
analytic. To introduce it we need some preliminaries. First of all, let us notice that 



( hplFl )-( hpW2 ) entail by simple maximum principle arguments (see [H Lemma 3.1]) 
that there exists a closed interval I\ C Iq such that any solution to (I2.25P satisfies 



Xoo E W 2 ' q (tt) Vg G [1, 00), Xoo(x) G h for all x E tt. 



Moreover, if we set, for v G V, 

E(v) :- 



-\Vv(x)\ 2 + W(v(x)) dx 



(2.26) 



(2.27) 



(which might be +00 if W(v) is not summable), there holds the following form of the 
Simon-Lojasiewicz inequality, which is a reformulation of [1, Prop. 4.4] (see also [21 
Prop. 4.2]): 
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Theorem 2.10. Let flhpWip - flEpTPgl ) hold and let 

W be real analytic on I . (2.28) 

Let Xqo be a solution to (12.251) . Then, there exist constants q, e > 0, ( G (0, 1/2), such 
that 

\E(v) - E^)] 1 ^ < c e \\Av + W'(v)\\v* (2.29) 

for all v G V such that 

\\v — Xco\\ VnC( j2) < e- (2.30) 

Remark 2.11. The above statement is given in a sligthly different fashion with respect 
to [H Prop. 4.4] since this version seems to be more suitable for our specific problem. 
However, we point out that our hypotheses entail those of [H Prop. 4.4]. Actually, by 
( 12.261) . it is clear that, as we possibly take a smaller e in condition fl2.30[) . then any v 
fulfilling (12.301) also satisfies [I, (4.5)]. Indeed, e can be taken so small that v ranges 
into Iq, where Lipschitz continuity of W holds (recall that Iq C I). 

Remark 2.12. Notice also that in this way we actually get rid of the possibly singular 
character of W at the boundary of I. Indeed, it is not excluded that there is a transient 
dynamics where W'(X(t)) may be unbounded. However, thanks to (I2.26P and the 
precompactness of the trajectory in C(fl), for sufficiently large times X(t) takes values 
into Iq, where W is bounded and analytic. As in [16], the key condition ensuring this 



property is ( hpW2 ), which essentially states that the leftmost and rightmost minima 
of W are interior to its domain. The gradient flow structure of the system entails that 
the solution eventually moves away from these minima. 

Here is, finally, our convergence result: 

Theorem 2.13. Let the hypotheses of Theorem 12. 71 hold. Furthermore, assume (12. 28ft 
and 

/■oo 

supt 1+s \\g(s)\\$.ds < oo, (2.31) 

t>0 Jt 

for some 5 > 0. Then, as (i?0j ^o) are initial data satisfying ( hp^p )-( hpX ), the uj-h'mit 
of the corrisponding trajectory (i?, X) of system (12. 101) -( [2. Ill) consists of a unique pair 
C#ooj ^oo), where is given by ( hpj'2 ) and Xoo G V is a solution to (12.251) . Moreover, 
as t y 00, ( 12.231) holds together with 

X(t)^Xoo weakly in H 2 (tt) and strongly in V DC (H) . (2.32) 

More precisely, if 

5 > — (2-33) 
1-2C 

where ( is as in (I2.29j) . then one can find t* > and a positive constant c* such that 

\X(t) -Xoo I < c*r^c, Vt>t*. (2.34) 
Otherwise, one can find Co G (0, C) so that 



(2.35) 



Colli - Hilhorst - Issard-Roch - Schimperna 



11 



a time t** > and a positive constant c„ such that 

\X(t) -Xoo| < c»t~^, \/t>t**. (2.36) 

Remark 2.14. Let us notice that (12.341) . (I2.36f) give the convergence rate of X(t) to 
with respect to the norm of H. Of course, an estimate of the rate of convergence 
in the norm of V could be obtained from the uniform bound corresponding to the last 
of (I2.17P and interpolation. 



3 A priori estimates and well posedness 

Let us first sketch an approximated version of system (I2.10l) -( j2.11l) along the lines of 
[29], [TU Sec. 3], to which we refer for more details. Namely, let us assume that j and 
W have been replaced in (I2.10l) -( l2.11l) by regularized functions j n and W n defined on 
the whole real line and such that 

jn, W n+2 "~* 3> W + 2 th<3 S6nSe ° f M ° SCO ^-l) 

(see, e.g., [3] for the definition of Mosco convergence and for the related notion of 
G- convergence of graphs). Moreover, we can assume [HI Sec. 3] that, for all n G N, 

W (r) 

W'(r) >-k Vr G R, — ^ > fi/2 Vr G R \ J , (3.2) 

r 

%(r) > a/2 Vr G R, (3.3) 



where /i, I , and a are as in (hpWl), ( hpVT2 ), ( hp j 1 ) (cf. [HI Sec. 3] for an example 



of a possible regularizing sequence {W„}). Then, we consider a family {(# n , X n )} ne ^ of 
(possibly local in time) solutions to the regularized problem specified by the subscript 
n. Existence of these solutions can be shown proceeding as in [291 Subsec. 3.1] and it 
turns out that (i? n , X n ), as well as u n := j' n ($ n ) an d W^(X n ), are regular enough to give 
a rigorous meaning to the forthcoming computations. Actually, we shall now deduce 
some a priori estimates with the aim of taking the limit as n /* oo. In this procedure, 
( 12.1 Op and (12. lip will be implicitly considered in their n-approximated form. Moreover, 
c will denote a positive constant, whose value is allowed to vary even inside one single 
row, but c may depend only on A, j, W, $o, X (and neither on n nor on t). When we 
need to fix the value of some specific c, we shall use the notation q, % > 0. The symbols 
C, Cj will denote constants which, instead, can explicitly depend on t G (0, oo), but do 
not explode as t \ 0. For simplicity, we shall proceed as if the solutions were defined 
for all times t > 0. Indeed, although this might be not true at the approximating 
level, it will certainly hold at the limit in view of the uniform in time character of the 
estimates and of standard extension arguments. 

Energy estimate. Test ( I2.10p by u n in the duality between V* and V and sum the 
result to (12.111) tested by X n ^ t in the duality of V* and V. The smoothness properties 
assumed on the approximating solutions, standard integration by parts formulas, and 
the cancellation of a couple of opposite terms then give 

^£„(0 n> Xn) + \Xn,t\ 2 + \\K\\v < \\\9\\v*, ( 3 - 4 ) 
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where we also used the Young inequality to split the duality product of g and u n 
resulting from the right hand side of (12.101) . Also, accordingly with (11.31) . we have set 

£„(#,*):= / {^\VX\ 2 + W n {X)+j n {d)). (3-5) 



Owing now to Qhp-$ p , ( hpXp ), and (I2.14p . we can integrate (13. 4p over (0, t) for arbitrary 



t > and deduce 



K||L2(0,t;V) + ||jn(^n)IU°°(0,t;Li(n)) < C, (3.6) 
\X n ,t\\L 2 (0,t;H) + ||^ n |U°°(0,t;V0 + || W n (X n ) \\ L^iO^L 1 (fi)) < C. (3.7) 



Let us note that, to obtain the second of (13. 7p . we used that, by (13.21) . 

3c, c > : J (^|Vw| 2 + W n (v)^j > c\\v\\ 2 v - c Vf 6 K 



(3.8) 



Second estimate. Let us now test ( 12.101) by $ n — $oo. Owing to the monotonicity 
of j' n and, more precisely, to (13. 3p . we can see that, both in the (IDirichletp and in the 
(IRobinP case, 

v*(Bu m $ n -$oo)v > 2c 1 ||tf n -tf 0O || 2 , (3.9) 

for some c\ > depending in particular on a. Thus, splitting the term depending on 
g by the Young inequality, we get 

\j t Wn ~ tfool 2 + CiR - VooWl < c\\g\\ 2 v « - v *{\\X n )X n)U $ n - i9 0O )v, (3.10) 
and the latter term is readily estimated as follows 

- v* (X\X n )X njt , $ n - ^oo) V = " (A'(X„)X n , t , d n - ^oo) 

< c(l + ||X n |U 4 (n))|^n,t| \\&n ~ #oo|U 4 (n) 

< c(l + \\X n \\ 2 v )\X n>t \ 2 + |||tf n -^oollv. ( 3 - n ) 

where we used ( |hpA[ ) and the continuous embeddings V, V C L 4 (Q). Hence, recalling 
( 13.71) and using ( 12.131) . we infer 

\\fln\\L~>(0,t;H) + \\$n ~ $oo |U 2 (0,t;V) < C - ( 3 -12) 

Moreover, arguing as in (13.111) . we obtain that the term \{X n ) t = A'(X„)X ni t in (I2.10p 
is uniformly controlled in L 2 (0,T;V*), so that (13.61) . (12.131) and the equality $ n)t = 
-X'(X n )X n)t - Bu n + g yield 

||tfn,*||l,2(0,T;V*) < C. (3.13) 

Third estimate. Thanks to ( |hpA| ), the continuous embeddings V, V C L A (Q), and 
(E3D (EIZD, we infer that 

||A / (X n )M„|| L 2 (0i4 . H) < c(l + ||Xn||L°°(0,t;L4(n)))||w n ||L2 (0it . X 4 (n)) < C. (3.14) 
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Then, we can test (12.1 ip by AX n and integrate over Q x (0, t). Using the first inequality 
of (13.21) . we have that 

(AX n , W' n {X n )) > -K\\X n f L ^ t . v) , (3.15) 

which can be controlled thanks to (13.71) . but only on bounded time intervals. Therefore, 
using also ( hpX ) it is not difficult to obtain 

||^n.|U 2 (0,t;H 2 (n)) + ||W / n(^n)IU 2 (0,t;H) < C, (3.16) 

where the second bound comes from a further comparison of terms in (12. lip . 

Limit asn-> +00 and existence. Standard compactness tools enable us to pass to 
the limit in the n-approximated versions of (I2.10p -( |2.1ip . Indeed, thanks to estimates 
(I3.6l) - (l3.7p . (I3.12p - (l3.13p . and (13.161) . there exist four limit functions 1?, u, X, v, defined 
from (0, 00) to H (at least), and a suitable subsequence of n (not relabeled) such that 
the corresponding subsequences {$ n }, {u n }, {X n }, {W^(X n )} fulfill 

$ n ^$ weakly star in iJ 1 (0, t; V*) R L°°(0; t; H) R L 2 (0, t; V), (3.17) 

u n -> u weakly in L 2 (0, £; V), (3.18) 

X n ^X weakly star in if 1 (0, t; H) fl L°°(0; t; V) R L 2 (0, t; H 2 (tt)), (3.19) 

-> v weakly in L 2 (0, t; H) (3.20) 

as n 00, for all t > 0. We note at once that the bounds in (I2.6p - (I2.9I) are certainly 
satisfied by the limit functions: in fact, it suffices to take the liminf in estimates (13.61) - 
fl3~7D . fl3~T2l - f[3~T3D which are uniform with respect to t. Next, fl3~T7]) . (EIinD, the Ascoli 
theorem and the Aubin compactness lemma (see, e.g., [321 Cor. 4, Sec. 8]) enable us 
to deduce that 

$ n -^$ strongly in C°([0; t];V*) R £ 2 (0, t; H), (3.21) 
X n ^X strongly in C°([0; t); H) R L 2 (0, t; V), (3.22) 

whence \'{X n ) — > A'(X) strongly in C°([0; £]; if) due to the Lipschitz continuity of A'. 
Then, (I3.18I) - (I3.19I) and the continuous embedding V C L 6 (f2) imply 

A'(X n )w n ->• A'(X)u weakly in L 2 (0, t; L 3/2 (fi)), (3.23) 
A(X n )t = V(X n )X„, t -»■ A(X) t = X'(X)X t weakly in L 2 (0, t; L 1 ^)). (3.24) 

At this point, we can pass to the limit in the n-approximated versions of (I2.10p -( |2.1ip 
to obtain (I2.10p (which a fortiori holds in V*) and 

X t + AX + v = X'(X)u inV*, (3.25) 

a.e. in (0, 00). Initial conditions (12.121) follow easily from (I3.21l) -( l3.22l) . Then, in order 
to conclude the existence proof, it remains to identify functions u and v, that is, to 
check that 

u = j'(ti), v = W'{X) a.e. in Vt x (0,oo). (3.26) 
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However, due to the Mosco convergences in (13.1 1) , it turns out that (cf., e.g., [3, 
Thm. 3.66]) the subdifferential operators f n , W' n + /tld G-converge to f, W + nld, 
as well as their extensions to L 2 (Q x (0, t)), for all t > 0. Then, we can apply the basic 
properties of G-convergence (see, e.g. [U Prop. 1.1, p. 42]) stating that if u n = j' n ($ n ) — ► 
u and d n — > $ weakly in L 2 (Q x (0,t)) and limsup ff nx ( 0t ) u n$n < ffnx(ot) u ^' then 
u = But, in our case this follows easily from (13.181) and (I3.2ip . The other iden- 

tification in (13.261) is a bit longer, since we first check that v + nX G (W f + Kld)(X) 
with the help of fl3~2"U|) and (l3~2"2l (recall that W' + kM is monotone by ( |hpP^l[ )), then 
extract the information v = W'(X). Hence, we conclude for the validity of ( 13.26ft . 

Uniqueness. Assume, by contradiction, that there exist two solutions ($j, Xj), i = 1,2, 
to problem fl2Tllil) -f l2T2l . let u t = f($i), i = 1,2, and set temporarily ■& := & x - d 2 , 
X := X\ — X2, u := u\ — u 2 . Now, we can take the difference of (12.101) written for the 
solutions corresponding to i = 1,2, integrate it in time from to t > 0, and test by 
u(t). At the same time, we test the difference of (12.1 ip by X(t) and sum the result 
to the previous relation. With the help of (hpjl) and (hpiyi) it is straightforward to 
deduce 

*\m\ 2 + ~ (\\f^(s)M\l + |X(*)| 2 ) + |VX(t)| 2 

<n\X(t)\ 2 + [ ((A'(X 1 )u 1 -A / (X 2 )M 2 )X-(A(X 1 )-A(X 2 ))u)(t). (3.27) 



To estimate the last term we use a remark from Kenmochi (see [IB] and, e.g. 
Lemma 3.2]): in fact, by the Taylor expansion and ( |hpA[ ) we have 

(A'(Xi)ui - \\X 2 )u 2 )X - (A(Xi) - \(X 2 ))u 

= u x {\{X 2 ) - A(Xi) - A'(Xi)(X 2 - Xi)) + u 2 (X(X 1 ) - A(X 2 ) - A'(X 2 )(X! - X 2 )) 
<^(\ui\ + \u 2 \)\X\ 2 . 

Therefore, the left hand side of ( 13. 27ft can be handled using the Holder inequality and 
the continuous embeddings V, V C L 4 (Q) in order to obtain 

a jT \$(s)\ 2 ds + \ (||/ ^(5)^|| 2 + |X(t)| 2 ) + jf \VX(s)\ 2 ds 

<c[ (l + \\ Ul {s)\\l +\\u 2 (s)\\ 2 v )\X(s)\ 2 ds + l- [ ||X(s)||^s. (3.28) 
Jo 1 Jo 

Finally, in view of the regularity ( 12. 7ft for u x , u 2 , the uniqueness property follows easily 
from the Gronwall lemma. This completes the proof of Theorem 12.21 

Proof of Theorem 12.61 Let us test ( 12.101) by u t in the duality between V* and 
V. Then, differentiate in time (12.111) . multiply the result by X t , and integrate over Q. 
Summing together the obtained relations, noting that a couple of terms cancel out, 
and using ( hpWl ), we infer 

Arlllnl" 
dt V2 



< -v*(9t,u) v + K\X t \ 2 + / X"(X)X 2 u. (3.29) 



n 
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Let us note that the computation above is just formal in the regularity setting of 
Theorem 12.21 However, the procedure might be made rigorous by working on the n- 
regularization sketched before and then passing to the limit. We omit the details just 
for brevity. By ( hp j 1 ) , we have 



{u U $t)>(T\M 2 - 



(3.30) 



Moreover, using QhpAD and once more the continuous embeddings V, V C L 4 (Q), we 
obtain 

1 



Finally, we note that 
Next, we set 



\"{X)Xiu < ^\\X t \\ l v + c\\u\\ z v \X t \ 



Jt,u) v < ||^|| v *||m||v- 



l.l 112 1, 



(3.31) 
(3.32) 
(3.33) 



and observe that, since g G L°°(0, oo; V*) thanks to Lemma [2.41 there exist constants 
C2, C3 > depending only on g and such that 

^\Hl + -\X t \ 2 <F+c 2 < \\u\\ 2 v + -\X t \ 2 + c 3 (3.34) 
for all t > 0. Thus, setting y := T + c 2 and using fl3~3UlHl3~m j32g) becomes 



+ a\$t\ 2 + ^|VX t | 2 < c|X t | 2 (l + \\u\\ 2 v ) + HIvlMlv 

<cy(i + y) + (i + \\u\\ 2 v){\\g t \\l 
< c (i + y)(i + y + \\g t \\ p v *). 

Then, noting that, since by the estimates fl2.7p - fl2.9p . 

ct+l 



1) 



sup 

t>0 Jt 



y(s) ds < 00, 



(3.35) 



(3.36) 



we can apply to 3^ the uniform Gronwall Lemma (see, e.g., [331 Lemma III. 1.1]), which 
yields ( 12. 16)) and the first bound in (12.171) . Next, integrating (13.351) in time over (i, t + 
1) for t greater than or equal to a given s > 0, we get the first of (12.151) and the 
second of (I2.17p . Observing that, by (hpjl )-(hpj2 ) (and the Poincare inequality in 
the ( IDirichletj) case), there exists c > such that 



IMIv > c\\d - ^oo||v, 
we get the second of (12.151) from (12.161) . Finally, by QhpAP we deduce 

ll^'(^)' U IU 00 (s,T;//) < c(l + ||X|| LO o( s T; y)) ||m|| L oo( s T;V) < C. 



(3.37) 



(3.38) 



Thus, using the first of (12.171) and (JhpWT| and viewing (12. lip as a time dependent 
family of elliptic equations with monotone (up to a linear perturbation) nonlinearities, 
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the last of (12.171) follows from the standard elliptic regularity theory. This, clearly, 
also gives (I2.18p . To conclude, we have to prove (12.201) in the case when (I2.19P holds. 
Coming back to the first row of ( 13.351) it then suffices to note that, for a.e. t > s, the 
right hand side is 

< c\X t \ 2 (l + \\u\\ 2 v ) + c|M|£ + c\\g t \\ q v , < c\X t \ 2 + c\\g t \\ q v * + c\\u\% , (3.39) 

where it is intended that, if q G (1, 2], then q* G [2, oo) denotes the conjugate exponent 
to g; otherwise, i.e., if q — 1, the latter term on the right hand side has to be omitted. 
Relation (I2.20p follows now by integrating (13.351) over (s, oo) and using (13.391) . Indeed, 
the terms on the right hand side of (13.391) are controlled, respectively, by the first of 
(12.91) . by (12.71) . (12.161) and interpolation, and by (12.191) . The proof of the Theorem is 
now complete. I 

4 Study of the u;-limit 

Proof of Theorem 12.71 Let us first notice that the convergence in (12. 24ft to some 
G H 2 (Q) is an immediate consequence of (I2.17P since X is weakly continuous from 

[s, oo) to H 2 (fl), whence the estimate ||X(t) 11^2^ < c holds true for all t G [s, oo). 
Analogously, the convergence $(t) — > in (I2.23P follows from (12.151) once one 

sees that flit) — > strongly in V* by ( 12. 6p . ( 12. 8p . and Lemma l2~4l indeed, due to the 

identification of the limit $00, it turns out that the entire family $(t) converges. 

Next, we show the convergence of u(t) in ( I2.23p . With this purpose, let us observe 

that, if fl2~2"2]) holds, then 

\\u t (t) \\ L i m < m) 1 < <i + n«(t) ii2*. (n) ) m) i 

<c(l+INlL-(l,ooiV))l^(*)l <c\#t(t)\- 

for a.e. t G (l,oo), say. Thus, by (12. 7p and (12.151) it is clear that 7 = u fulfills the 
assumptions of Lemma [2.41 with p = 2 and X = L 1 (fi) and consequently u(t) — > in 
L x (f2) as t y 00. Then, to show that (12.231) holds (i.e., u(t) weakly converges in V), it 
is now enough to point out the bound in ( 12.161) . 

Let us consider, instead, the case when ( 12.211) holds. Under this condition, using 
(|2~7D . (jZHD and Remark we modify fl3~3D as 

< - A \\u\\ 2 v + i|X t | 2 < Z := T + \\gf v , < \\u\\ 2 v + ^\X t \ 2 + ^\\g\\ 2 v *. (4.1) 

Thus, in view of the first line of ( 13.351) . it is not difficult to deduce 

^-Z < c(Z 2 + 1 + \\g\\ 2 v , + \\g t \\ 2 v *) a.e. in (0,oo), (4.2) 

whence the convergence u(t) — > with respect to the strong topology on V is a con- 
sequence of [311 Lemma 6.2.1, p. 225]. Then, recalling (13.371) . in this case it happens 
that (12.231) is improved into 

u(t) — > 0, $(£)—>• strongly in V. 
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To conclude, it remains to show that the limit value X^ in (I2.24p satisfies §ZZ 
and this can be done in a completely standard way. Namely, defining (X n , u n ) : (0, 1) — ■> 
V x V as (*„,«„)(■) := (X,«)(t„ + •)> from (EH]) and (ETJ) it is clear that 

(X n ,«n) ->(*«,, 0) strongly in C°([0, 1]; if) x L 2 (0, 1; V). (4.3) 

as n y oo. Here, of course we also used the strong convergence of X n<t to in 
L 2 (0, 1; if). Then, by passing to the limit in 

X ntt + AX n + W'(X n ) = X'(X n )u n inV*, a.e. in (0,1), (4.4) 

it is not difficult to check that Xoo solves the stationary problem. Indeed, owing to 
(I2.17p . weak star compactness, and (14.31) it turns out that X n — + X^ weakly star in 
f°°(0, 1;H 2 (Q)), whence AX n -> AX^ weakly star in L°°(0, l;if). Moreover, using 
(12.1 8p and exploiting the maximal monotonicity of W + k Id (which is a continuous 
and increasing function thanks to (hpWl)), one verifies that W'(X n ) tends to W\X 00 ) 
weakly star in L°°(0, l;if) (and then weakly in L 2 (Q x (0,1))) with the help of the 
strong convergence X n —>■ X^ in L 2 (Q x (0,1)) and of [H Prop. 1.1, p. 42]. Finally, 
thanks to (14. 3p . QhpAP and the continuous embedding V C L 6 (Q), we infer that 

\'{X n )u n -> strongly in f 2 (0, 1; L 3/2 (Q)), 

and consequently the right hand side of (14.411 tends to in f 2 (0, 1; V*). This concludes 
the proof of Theorem 12 .71 I 

Proof of Theorem 12.131 We proceed partly as in [7J Sec. 3], [T7]. Let us assume, 
for simplicity, that ( 12.331) holds (otherwise, we can replace ( with a value Co such that 
( 12.3511 is satisfied and notice that Theorem 12.101 still holds with (o in place of £ in 
( 12.2911 ) . Letting Xoo be an element of the cu-limit of X(-), we can set (cf. Theorem 12.101 
for the notation) 

S := {t > : \\X(t) - Xoo||ynco(n) < e/3}. (4.5) 
Clearly, £ is unbounded. Next, for t e S, we put 

r(t) := sup {t' > t : sup \\X(s) - ^oo||y n c°(fI) < e } ( 4 - 6 ) 
se[t,f] 

and observe that, by continuity, r(t) > t for all t E S. Let us fix t E S and divide 
J7 := [to? r (^o)) (where r(t ) might well be +oo) into two subsets: 

A x := {t E J : \X t (t)\ + |K*)||y > ( J \\g(s)\\ 2 v * ds) j, (4.7) 
A 2 :=J\A 1 . (4.8) 

Next, we define (cf. (l2~27D ) 

$(t) := J JWx,t)) dx + ~ ^" T(t0) ||^( S )|| 2 « ds + £(*(*)) - S(Xoc). (4.9) 
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Then, it is not difficult to see that 

* , (*)<-(^ll«(*)llv + l^(*)l a )- ( 4 - 10 ) 

We remark that $ is absolutely continuous thanks to fl2.15p - fl2.18p and [3J, Lemme 3.3, 
p. 73]. This justifies the above computation for a.e. t G J . We then have (cf. also [211 
(3-2)]) 

^(^Nign^W^-Cl^l^^lkWII^ + W)! 2 ). (4.11) 

Now, let us estimate $ from above. If t G A%, thanks to ( 12. 27ft and Theorem 12.101 we 
obtain 



m)\ 1 ~ i <\E{x{t))-E{x QO )\^+ j(m 



i-C 



\9W\\v 



2 ds 



i-C 



i-C 



<c t \\ -X t (t)+\'(X(t))u\\ v * + +l^(*)l + IK*)l|v, (4.12) 



where we also used ( 14.71) . Note now that, by ( ]hpA[ ), the last of (I2.17p . and well-known 
continuous embeddings, we have 

c t \\ - X t (t) + \'(X(t))u\\ v * < c(\X t (t)\ + \\u(t)\\ v ). (4.13) 

Moreover, by convexity of j, j(i?oo) — 0, and Holder's inequality, we infer 

0< [ j{ti{x,t)) dx < [ j'(tf(x,t))(tf(x,t) -^oo) dx < \u(t)\\&{t) -i?oo|- (4.14) 



Using once more (I3.37P together with (12.161) and recalling that ( G (0, 1/2), one gets 



Km 



i-C 



< c\\u 



2(1-0 
V 



< c\\u 



Thus, collecting (j02jl - (O51) . we finally have 

\m\ l ~ C <c(\Xt{t)\ + \\u{t)\\v), 
whence from (14.111) we obtain 

ll^)llv + |^(t)|<-^^(|$Nign$)(t). 



(4.15) 



(4.16) 



(4-17) 



Since $ is decreasing by ( I4.10p . integration in time entails that ||lt||y and \X t \ are 
summable over A\. Of course, the same holds over A 2 by (12.311) and (14. 7ft - (14 .8 1) . Thus, 
we conclude that X t G L 1 (j7;i?). 

From this point on, the proof proceeds exactly as in [T71 Sec. 3] . Namely, a simple 
contradiction argument yields that r(i ) = oo as t G S is sufficiently large. This 
implies that X t G ^(to, +oo; H), whence the convergence (in H) of the whole trajectory 
X(t) to Xqo follows. More precisely, this convergence holds strongly in V fl C(O) by 
precompactness of the trajectory (cf. ( 12. 17ft ). Finally, the technical argument of [T71 
Sec. 3] leading to estimate (I2.34p (or (12.361) ) can be repeated just by adapting the 
notation. I 



Colli - Hilhorst - Issard-Roch - Schimperna 



19 



References 

[1] S. Aizicovici and E. Feireisl, Long-time stabilization of solutions to a phase-field 
model with memory, J. Evol. Equ., 1 (2001), 69-84. 

[2] S. Aizicovici, E. Feireisl, and F. Issard-Roch, Long time convergence of solutions 
to a phase-Geld system, Math. Methods Appl. Sci., 24 (2001), 277-287. 

[3] H. Attouch, "Variational Convergence for Functions and Operators" , Pitman, Lon- 
don, 1984. 

[4] V. Barbu, "Nonlinear Semigroups and Differential Equations in Banach Spaces", 
Noordhoff, Leyden, 1976. 

[5] H. Brezis, "Operateurs Maximaux Monotones et Semi-groupes de Contractions 
dans les Espaces de Hilbert", North-Holland Math. Studies 5, North-Holland, 
Amsterdam, 1973. 

[6] G. Caginalp, An analysis of a phase field model of a free boundary, Arch. Rational 
Mech. Anal., 92 (1986), 205-245. 

[7] R. Chill and M.A. Jendoubi, Convergence to steady states in asymptotically au- 
tonomous semilinear evolution equations, Nonlinear Anal., 53 (2003), 1017-1039. 

[8] P. Colli, G. Gilardi, E. Rocca, and G. Schimperna, On a Penrose-Fife phase-field 
model with nonhomogeneous Neumann boundary conditions for the temperature, 
Differential Integral Equations, 17 (2004), 511-534. 

[9] P. Colli, M. Grasselli, and A. Ito, On a parabolic-hyperbolic Penrose-Fife phase- 
Geld system, Electron. J. Differential Equations, 2002, No. 100, 30 pp. (electronic) 
[Erratum, Electron. J. Differential Equations, 2002, No. 100 erratum, 32 pp. 
(electronic)]. 

[10] P. Colli and Ph. Laurengot, Weak solutions to the Penrose-Fife phase Geld model 
for a class of admissible heat Gux laws, Phys. D, 111 (1998), 311-334. 

[11] P. Colli, Ph. Laurengot, and J. Sprekels, Global solution to the Penrose-Fife phase 
Geld model with special heat Gux law, in "Variation of domains and free-boundary 
problems in solid mechanics", P. Argoul, M. Fremond & Q.S. Nguyen eds., Solid 
Mech. Appl., 66, Kluwer Acad. Publ, Dordrecht, 1999, pp. 181-188. 

[12] P. Colli and P. Plotnikov, Global solution to a quasistationary Penrose-Fife model, 
Indiana Univ. Math. J., 54 (2005), 349-382. 

[13] E. Feireisl, F. Issard-Roch, and H. Petzeltova, Long-time behaviour and conver- 
gence towards equilibria for a conserved phase Geld model. Partial differential 
equations and applications. Discrete Contin. Dyn. Syst., 10 (2004), 239-252. 

[14] E. Feireisl and G. Schimperna, Large time behaviour of solutions to Penrose-Fife 
phase change models, Math. Methods Appl. Sci., 28 (2005), 2117-2132. 



20 LONG TIME CONVERGENCE FOR PHASE FIELD MODELS 

[15] E. Feireisl and F. Simondon, Convergence for semilinear degenerate parabolic 
equations in several space dimensions, J. Dynam. Differential Equations, 12 
(2000), 647-673. 

[16] M. Grasselli, H. Petzeltova, and G. Schimperna, Long time behavior of solutions 
to the Caginalp system with singular potential, Z. Anal. Anwendungen, 25 (2006), 
51-72. 

[17] M. Grasselli, H. Petzeltova, and G. Schimperna, Convergence to stationary solu- 
tions for a parabolic-hyperbolic phase-field system, Comm. Pure Appl. Anal., 5 
(2006), 827-838. 

[18] N. Kenmochi, Uniqueness of the solution to a nonlinear system arising in phase 
transition, in "Nonlinear analysis and appplications," N. Kenmochi, M. Niezgodka, 
and P. Strzelecki eds., GAKUTO Internat. Ser. Math. Sci. Appl., 7, Gakkotosho, 
Tokyo, 1995, pp. 261-271. 

[19] A. Ito and N. Kenmochi, Inertial set for a phase transition model of Penrose-Fife 
type, Adv. Math. Sci. Appl., 10 (2000), 353-374 [Erratum, Adv. Math. Sci. Appl., 
11 (2001), 481]. 

[20] A. Ito, N. Kenmochi, and M. Kubo, Non-isothermal phase transition models with 
Neumann boundary conditions, Nonlinear Anal., 53 (2003), 977-996. 

[21] M.A. Jendoubi, Convergence of global and bounded solutions of the wave equation 
with linear dissipation and analytic nonlinearity, J. Differential Equations, 144 
(1998), 302-312. 

[22] N. Kenmochi and M. Kubo, Weak solutions of nonlinear systems for non- 
isothermal phase transitions, Adv. Math. Sci. Appl., 9 (1999), 499-521. 

[23] S. Lojasiewicz, Une propriete topologique des sous ensembles analytiques reels, 
Colloques du CNRS, Les equations aux derivees partielles, vol. 117, 1963. 

[24] O. Penrose and P. C. Fife, Thermodynamically consistent models of phase-field 
type for the kinetics of phase transitions, Phys. D, 43 (1990), 44-62. 

[25] O. Penrose and P. C. Fife, On the relation between the standard phase- field model 
and a "thermodynamically consistent" phase-field model, Phys. D, 69 (1993), 107- 
113. 

[26] E. Rocca and G. Schimperna, Universal attractor for some singular phase transi- 
tion systems, Phys. D, 192 (2004), 279-307. 

[27] E. Rocca and G. Schimperna, Universal attractor for a Penrose-Fife system with 
special heat flux law, Mediterr. J. Math., 1 (2004), 109-121. 

[28] E. Rocca and G. Schimperna, Global attractor for a parabolic-hyperbolic Penrose- 
Fife phase field system, Discrete Contin. Dyn. Syst., 15 (2006), 1193-1214. 



Colli - Hilhorst - Issard-Roch - Schimperna 



21 



[29] R. Rossi, Well-posedness and asymptotic analysis for a Penrose-Fife type phase 
Geld system, Math. Methods Appl. Sci., 27 (2004), 1411-1445. 

[30] W. Shen and S. Zheng, Maximal attractors for the phase-field equations of Penrose- 
Fife type, Appl. Math. Lett., 15 (2002), 1019-1023. 

[31] L. Simon, Asymptotics for a class of non-linear evolution equations, with applica- 
tions to geometric problems, Ann. Math., 118 (1983), 525-571. 

[32] J. Simon, Compact sets in the space L p (0,T;B), Ann. Mat. Pura Appl. (4), 146 
(1987), 65-96. 

[33] R. Temam, "Infinite Dimensional Dynamical Systems in Mechanics and Physics" , 
Springer- Verlag, New York, 1988. 



[34] S. Zheng, "Nonlinear Evolution Equations", Chapman Hall/CRC Monogr. Surv. 
Pure Appl. Math., 133, Chapman k Hall/CRC, Boca Raton, FL, 2004. 



